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Abstract 

Mass number and isospin dependence of symmetry energy coefficients for isobaric chains with 
A = 8, 12, 16, 20, • • • , 100 have been globally studied in relativistic mean-field (RMF) theory with 
the PK1 effective interaction and compared with the corresponding pseudo data from the experi- 
mental binding energy and semi-empirical mass formula. The mass dependence of symmetry energy 
coefficient both from the RMF calculation and pseudo data for different isospin values has been 
fitted with and without the surface term b v j A^l^ . The results confirm that only if with the surface 
term in the fitting of symmetry energy coefficient in finite nuclei can one get the consistent value 
as that in infinite nuclear matter. 
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I. INTRODUCTION 



In the past decades, many exotic nuclei far from the valley of /3-stability line with extreme 
N/Z ratio have been produced with the new generation of radioactive- ion beam facilities 
These exotic nuclei provide us useful information on the equation of state (EOS) of asymmet- 
ric nuclear matter, which has not been well determined and is important for understanding 
both the structure of unstable nuclei and the properties of neutron stars. One of the most 
important quantities is the nuclear symmetry energy, which affects significantly the binding 



energy and rms radii of neutron-rich nuclei 



2|. Furthermore, the chemical composition, 



the evolution of lepton profiles, cooling process and the neutrino fluxes in neutron stars 
depend strongly on the nuclear symmetry energy 4]. Therefore, more and more attention 
has been payed to study the symmetry energy in both infinite nuclear matter and finite 



nuclei 15 



6. 
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In infinite asymmetric nuclear matter, the energy density is usually expanded in terms of 
isospin asymmetry (3 = (p n — P P )/p as £(p,/3) = S(p,/3 = 0) +a' sym (p) (3 2 +0(/3 2 ), where the 
nuclear symmetry energy coefficient is given by a' sym (p) = — — 7^ — ^=0, and £(p,f3 = 0) 
is the energy density in symmetric nuclear matter. Many investigations have already been 
carried out to study the density and isospin dependence of symmetry energy coefficient 
a' sym Q. 

In finite nuclei, nuclear symmetry energy has been studied in the finite-range droplet 
model (FRDM) [5J, Random Phase Approximation (RPA) based on the Hartree-Bogoliubov 



approach 



13|, and energy density functionals of Skyrme force (9] as well as relativistic 



meson-exchange force IOJ. It has been found that the nuclear symmetry energy coefficient 
a sym extracted from E sym = a sjm T(T + 1) is in better agreement with the empirical data 
than that from E sym = a sym T 2 . 

In this paper, the mass number and isospin dependence of symmetry energy coefficients for 
isobaric chains with A = 8, 12, 16, 20, • ■ • , 100 will be globally studied in a covariant density 
function theory and compared with the corresponding pseudo data from the experimental 
binding energy and semi-empirical mass formula. The paper is organized as follows. In 
Sec. [Til we briefly describe the method used to obtain the nuclear symmetry coefficient in 
finite nuclei. The results and discussion will be given in Sec. IIHI Finally, a short summary 
is presented in Sec. [IVl 
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II. THE METHOD 



The conventional semi-empirical mass formula of Bethe and Weizsacker in liquid droplet 
model (LDM) for the binding energy of a nucleus is as follows fl^ |. 

Z 2 



B{N,Z) = a v A-a s A 2/3 



a. 



A 1 / 3 



(N-Z) 2 5 m 

-a-i ^ + a PAt/i' W 

with N and Z the number of neutrons and protons, A the total number of nucleons, A = 
N + Z. a v , a s , a c , and aj are positive and represent the coefficients of bulk or volume, surface, 

Coulomb, and symmetry energies respectively, and a p is a simple parametrization of pairing 

(-1)*+ (-l) z 
with 5 = ± — J —-± — '—. 
2 

The research of interest in this work is the symmetry energy term in Eq.flT]), which can 
be rewritten in terms of the isospin value T, 

E sym (A, T) = a 7 (iV ~ Z)2 = a sym (A, T)T 2 , (2) 

where t 



wnere the isospin value of a nucleus is given by T = \T Z \ = \N — Z\/2. However, recent 
studies , [lol . [13! . I5I ] showed that it is more proper to parameterize the nuclear symmetry 
energy into following form, 



E sym (A, T) = a sym (A, T)T(T + 1), (3) 

which will be adopted to extract the symmetry energy coefficient subsequently 

On the other hand, according to the LDM formula (CQ), one notices that the difference in 
binding energies of isobaric nuclei with same odd-even parity is only related to the Coulomb 
energy and symmetry energy terms, 

B(A,T) - B(A,T = 0) 
= [E CouL (A,T) - E Coul .{A,T = 0)} 

+ [E sym (A,T) - E sym (A, T = 0)]. (4) 

Therefore, the pseudo experimental value of symmetry energy in Eq.Q is determined by 
subtracting the Coulomb energy from the experimental data of nuclear total binding energy 
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c.lS . 



E sym (A, T) 

= [B(A,T) - B(A,T = 0)] - [Ecoui.(AT) 

-E Coul .(A,T = 0)}, (5) 

where the Coulomb energy E CouL (A, T) of a specific nucleus is calculated by assuming the 

nucleus as an uniformly charged spherical droplet, i.e., 

, m\ 3Z¥ r 0.7636, 
Ecoui.(AT) = - — [1 - ^r]. (6) 

The charge radius is given by R c = 1.2 A 1 / 3 . The second term is from the Fock-exchange term 
of proton. Even if there is quadrupole deformation (3 (typically about 0.1 — 0.5) deviating 
from sphere in a nucleus, it just reduces the Coulomb energy by a factor of /3 2 /4ir, which is 
negligible. Therefore, it is justified to consider the nucleus as a spherical liquid drop rather 
than a deformed one. Finally, one obtains the pseudo experimental value of symmetry energy 
coefficient from Eqs.(j3]) and (jSJ), which depends on the mass number A and isospin value T. 

The theoretical study of symmetry energy coefficients for finite nuclei is based on a 
covariant density functional theory with nucleon-nucleon interacting by exchange of various 
mesons and photo for nuclear ground-state energies. The non-linear version of relativistic 
mean-field theory is used. The corresponding effective Lagrangian density are composed of 
the free parts of nucleon field (ip), two isoscalar meson fields (a and u^), isovector meson 
field (Pfj), photon field (A^) as well as their couplings [l6,[l7, 18], 



C = ip [iYdf, -m- g a a - g^uj^ - g^r ■ 
1 



2 -e(l - r 3 )7% 



-fir ■ iv + u p fr) - \f^f^ (7) 



where m and rrii(gi) (i = a, u^, p M ) are masses (coupling constants) of the nucleon and the 
mesons respectively and 

QT = &>u v - d v u", (8a) 
fir = d 1 *? - d"p», (8b) 
pw> = d»A v -d u A» (8c) 



are field tensors of the vector mesons and the photon field. We adopt the arrows to indicate 
vectors in isospin space and bold types for the space vectors. Greek indices \i and v run 
over 0, 1, 2, 3, while Roman indices i, j, etc. denote the spatial components. 

The nonlinear self-coupling terms U a (o), U^u^), and U p (p^) for the cr-meson, cu-meson, 
and p-meson in the Lagrangian density (0) respectively have the following forms: 

U a (<r) = ^rn 2 y + -g 2 a 3 + -g 3 a 4 , (9a) 
Uuiwn) = -mlu^Uft + ^c 3 {u^u^ 2 , (9b) 
U P (P») = Imlp^Pv (9c) 

From Eq.([71) combined with the mean- field approximation, one obtains the corresponding 
energy density functional -Ermf = (^l-^l^)) where H is the corresponding hamiltonian 
operator. The equations of motion for the nucleon and the mesons can be obtained by 
requiring that the energy functional be stationary with respect to the variations of single- 
particle density p and meson fields. Finally, one obtains a Dirac equation for nucleon as well 
as inhomogeneous Klein-Gordon equations for meson fields and photon field. The source 
terms for Klein-Gordon equations are sums of bilinear products of Dirac spinors restricted 
in the Fermi sea. 

With the restriction of spherical symmetry, both the Dirac equation for the nucleon and 
the Klein-Gordon equations for the mesons and the photon become radially dependent only. 
These coupled equations are solved iteratively using the shooting method combined with 
Runge-Kutta algorithms and employing appropriate boundary conditions in a spherical box 
of radius R = 20 fm with a step size of 0.1 fm. More details about the numerical techniques 
can be found in Ref. [l8|. 

Moreover, the proper treatment of center of mass (cm.) motion has been found very 
important in the explanation of Coulomb energy displacement in mirror nuclei [jj3]. We 
adopt the same cm. correction to the total energy after variation, as it has been used in 
adjusting the parameter set PK1 j^ . 

£ ™» = "aSA<^>- < 10 > 

where m is the mass of neutron or proton, P cm = Pi is the total momentum in the cm. 
frame. The total energy for the nuclear ground state is given by the summation of nuclear 
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part -Ermf and cm. correction, 



Etot. = E RM f + E™ c . (11) 

The nuclear symmetry energy for certain A and T is finally obtained as follows, 

E^ e m (A,T) = E tot .(A,T = 0) - E tot .(A,T) (12) 

where E to t.(A,T = 0) is the total energy of N = Z isobaric nucleus with mass number A. 
The theoretical value for the coefficient of symmetry energy a sym (A,T) is thus determined 
hyE^ m (A,T)/T(T+l). 



III. RESULTS AND DISCUSSION 



In this section, we perform systematical calculations for A = 8, 12, 16, 20, • ■ ■ , 100 isobars 
with the effective interaction PK1. The main interesting is focused on the coefficient of 
symmetry energy, which is only associated with symmetry energy. It is justified to switch off 
the Coulomb interaction in the particle-hole channel of nucleon-nucleon effective interaction 
in the effective Lagrangian density flTj). 
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FIG. 1: The coefficients of symmetry energy for nuclei with T = 2,4, • • • , 12 as a function of the 
corresponding mass number A. The filled symbols represent the pseudo experimental data and the 
open symbols are the RMF predictions. 
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Figure [T] shows the predicted nuclear symmetry energy coefficients a sym in nuclei with 
T = 2, 4, ■ • • , 12 as a function of the corresponding mass number A in comparison with 
the corresponding pseudo data. It shows that good agreement in a sym has been achieved 
between RMF calculation and the pseudo data except the very light nuclei with T = 2. The 
observed strong oscillating in a sym of RMF calculation for small T are associated with shell 
closures |lo| . In contrary, there is no effect of shell structure in the LDM formula, which 
gives a moderate changing of function of mass number A. Furthermore, for larger 

value of T, the symmetry energy becomes less sensitive to the shell structure and its value 
is stabilized. As the consequence, a better agreement between RMF calculation and the 
pseudo data is found for large isospin value T. 

In the following, we parameterize the coefficient of symmetry energy into following two 
forms: 

A A 4 / 3 

where a' syra in the first form corresponds to the standard definition of coefficient of symmetry 
energy in infinite nuclear matter, while b v and b s are the so-called volume and surface com- 
ponents in symmetry energy of finite nuclei and its physical origin is traditionally explained 



in terms of the kinetic energy and mean isovector potential contributions [2l|. Duflo and 
Zuker j^J performed a global fit to nuclear masses and found b v = 134.4, b s = 203.6. Ban et 
al. [loj], made a least-square fit to the calculated nuclear symmetry energy for the A = 40, 
48, 56, 88, 100, 120, 140, 160, 164, and 180 isobaric chains and found b v = 133.2, b s = 220.3. 

We carry out a global fit of the obtained coefficients of symmetry energy from pseudo 
data and RMF calculation via the Levenberg-Marquardt method. Table [J presents the 
parameter a' sym in coefficients of symmetry energy a sym for isospin values of T = 2, 4, 6, 
and 8 respectively. It shows that the parameters a' sym from both the pseudo data and RMF 
calculation are more than two times smaller than the infinite symmetric nuclear matter value 
37.6 MeV for PK1. 

Table [Til presents the parameters b v and b s in coefficients of symmetry energy for isospin 
values of T = 2, 4, 6, and 8 respectively. It shows that the parameters b v and b s are consistent 
with those given in Refs. 0, [loj and the coefficient of symmetry energy in infinite nuclear 
matter a' syrn = b v /4 with A going to infinity in Eq.( fT3l . It confirms the conclusion that only 



TABLE I: The parameter a' sym = 4:a' sym /A in coefficients of symmetry energy a sym for isospin values 
of T = 2, 4, 6, and 8 respectively. 



r 


2 


4 


6 


8 


Exp. 


13.5 


16.8 


18.3 


19.2 


Theory 


9.2 


14.6 


16.0 


18.2 



TABLE II: The parameters b v and b s in coefficients of symmetry energy a sym — b v I A - b s /A^ 3 
for isospin values of T = 2, 4, 6, and 8 respectively. 





Exp. 


Theory 


T 


b v 


b s 


b v 


b s 


2 


131.4 


183.9 


134.6 


232.5 


4 


92.2 


68.4 


111.3 


157.8 


6 


121.5 


175.9 


137.0 


244.5 


8 


128.1 


206.1 


138.9 


254.8 



with the surface term b v /A 4 ^ 3 in the fitting of symmetry energy coefficient can one get the 
consistent value as that in infinite nuclear matter. 

Figure [2] shows the coefficients of symmetry energy for nuclei with A = 40, 60, 80 and 100 
as a function of the corresponding isopsin value T in comparison with the corresponding 
pseudo data. It shows that the coefficients of symmetry energy a sym for fixed mass number 
decrease gradually with the isospin value T. 

IV. SUMMARY 

The mass number and isospin dependence of symmetry energy coefficient in finite nuclei 
have been studied in RMF theory with the effective interaction PK1 and compared with 
the pseudo data from the experimental binding energy and semi-empirical mass formula. 
It has been found that the nuclear symmetry energy coefficient obtained from the RMF 
calculations decreases gradually with the mass number and isospin value, which is in good 
agreement with the pseudo data. The mass dependence of symmetry energy coefficient 
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FIG. 2: The coefficients of symmetry energy for nuclei with A = 40, 60, 80 and 100 as a function 
of the corresponding isopsin value T = (N — Z)/2. The filled symbols represent the pseudo 
experimental data and the open symbols are the RMF predictions. 

for different isospin values has been fitted with and without the surface term b v /A i ^ 3 . The 
results confirm that only if with the surface term in the fitting of symmetry energy coefficient 
in finite nuclei can one get the consistent value as that in infinite nuclear matter. 
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